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1. Introduction 



This is a report on joint work with Alan Rendall |AR.QC| ], on non-oscillatory 



singularities in four dimensional space-times with scalar field or stiff fluid matter. 



Before considering the results of [AROC] in detail, I will discuss the BKL proposal 
which motivated the work. 

The singularity theorems of Penrose and Hawking guarantee the existence of 
spacetime singularities under reasonable assumptions, but give little information 
about the singularities they predict. The most detailed proposal for the structure 
of space-time s i ngularit ies is due to Belinskii, Khalatnikov and Lifshitz (BKL), see 



|LK63 BKL7C , BKL82] and references therein. According to the BKL proposal, 
singularities in generic four dimensional space-times are space-like and oscillatory, 
while generic space-times with stiff fluid, including massless scalar fields, have sin- 



gularities which, according to a proposal of Belinskii and Khalatnikov, see | BK72 |, 
are space-like and non-oscillatory. Further, generic space-times will, according to 
this picture, have the property that spatial points decouple near the singularity, and 
the local behavior is asymptotically like spatially homogeneous (Bianchi) models 
locally near the singularity. In particular, in the non-oscillatory case, spatial deriva- 
tives become unimportant, the term "asymptotically velocity term dominated" is 
often used to describe this phenomenon. Space-times with non-stiff matter appear 
to have the property that asymptotically close to the singularity, matter becomes 
insignificant and the main features of the evolution is determined by curvature. On 
the other hand, in the case of stiff matter, the picture is that the matter becomes 
dominant near the singularity, leading to non-oscillatory behavior. We will refer to 
the collection of ideas described above as the BKL proposal. 

The Bianchi models arc space-times with spatially locally homogeneous geom- 
etry. This is an interesting toy model for the structure of singularities, and the 
so-called mixmaster universe, Bianchi IX, has been used as a model for the typical 
local behavior, near the singularity, of gravity in four dimensions, with non-stiff 
matter. In the case of Bianchi space-times, the Einstein equations are a system 
of ODE's and the asymptotic behavior of the space-time at the singularity can 
be studied using dy namical systems methods. It has recently been proved by H. 



Ringstrom [RinOOb] that the singularity in generic Bianchi class A space-times 
with non-stiff matter, which includes Bianchi I X, is os cillatory. In the stiff fluid or 
scalar field case, Bianchi IX is non-oscillatory | Rin004 §20]. In either case, curva- 



ture blows up at the singularity, which has as a consequence that the space-time 
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is inextendible, i.e. the conclusion of the cosmic censorship conjecture cf. |Pen82|, 
holds for Bianchi class A. 

For generic space-times, in the case of non-oscillatory singularities, the BKL 
proposal leads to the space-time being asymptotically Kasner at the singularity, 
with the Kasner exponents depending on the "point on the singularity". In the 
oscillatory case, the proposal indicates that locally near the singularity, the space- 
time returns infinitely often to a state which is approximately a non-flat Kasner. 
From this it would follow that certain curvature scalars such as the Kretschmann 
scalar k = ^ Rapyg^ R a ^ jS become unbounded near the singularity. Thus, verify- 
ing the correctness of the BKL proposal would go a long way towards solving the 
cosmic censorship problem. 

Barrow [Bar78| discussed the structure of the early universe for space-times with 
stiff matter and used the term "quiescent cosmology" for this picture. Here we will 
use the term "quiescent singularity" for a singularity with non-oscillatory behavior 
induced by the pre sence o f matter. In section |^ below I will report on joint work 
with Alan Rendall | AROO | , where we construct families of four dimensional space- 
times with stiff fluid or scalar field matter, with quiescent singularities. These 
families have the same number of free functions as the full space of solutions to the 
Einstein equations, and therefore this result supports the BKL proposal described 
above for this case. 

There are important situations where classes of space-times with symmetry ex- 
hibit non-oscillatory behavior at the singularity, even without the presence of stiff 
matter. An example is the Gowdy class of space-times where both numerical and 
analytical work gives support to the notion that generic Gowdy space-times have 
non-oscillatory singularities. Gowdy space-times have a 2-dimensional symmetry 
group G2, with the action of G2 generated by spatial Killing fields with vanishing 
twist. Rendall and Kichcnassamy used Fuchsian methods to construct families, 
with the "right" number of degrees of freedom, of Gowdy space-times with non- 
oscillatory singularity. This result uses a singular version of the Cauchy-Kowalewski 
theorem, and thus requires real analyticity for the data. The Fuchsian method is 
also used in the four dimensional results described in section |2|. The analyticity 
condition has recently been removed for the Gowdy case by Rendall jRcnOOH . In 
general, G2 space-times with matter or with non-vanishing twist have oscillatory 
singularities ]WIB98 |. 

The dimension of the space-time is significant. It has been argued, using a BKL 
type analysis, that in space-ti me dime nsion D > 11, vacuum gravity generically 
has non-oscillatory singularity [ DHS85 1 . A low energy limit of string theory leads 
to consideration of Einstein equations in D dimensions, D = 10, 11, with matter 
consisting of axions, dilatons (scalar field) and form fields. This model is some- 
times called "stringy gravity" . The possibility that these space-times have non- 
oscillatory singularities was exploited in the so-called pre-big-bang scenario, where 
one views the present universe as arising out of a (future) cosmological singularity 
in a previous universe, the present universe being blown up to a scale determined 
by the dilaton. This was studied by Buonnanno et al ]BDV99 | where it was shown 
using formal expansions that the singularities of stringy gravity in spherical symme- 
try are non-oscillatory. Families of stringy gravity solutions with Gowdy symmetry 
and non-oscillatory singularity have been constructed using the Fuchsian method 
|NTM0C ] . Recent work by Damour and Henneaux |DH00b, DHOOa] indicates that 
stringy gravity without extra symmetry does in fact have oscillatory singularities. 



QUIESCENT COSMOLOGICAL SINGULARITIES 



3 



2. Quiescent singularities 

Here we describe the main results of [ AR0(| . For simplicity we restrict the discus- 
sion here to the case of massless scalar field matter. We consider four-dimensional 
space-times with metric of the form 

ds 2 = -dt 2 + g ab (t)9 a <8> 9\ 

so that tisa Gaussian time coordinate and {0 a } is a coframe dual to a frame {e a }. 
In this case, the second fundamental form of a hyper-surface t ^constant is given 
by k ab = —\dt9ab- The Einstein field equations are 

^R Qfi = 8nV a (t>V <t>, 
which can be written in the following equivalent 3 + 1 form. The constraints are: 
R - k ab k ab + (trfc) 2 = 8ir[(d t <f>) 2 + V a 0V a 0], 
V a /c afc - Vfctrfc = -8Trd t <f)V b <f), 
and the evolution equations are: 

dtgab = -2k a b, 

d t k a b = R a b + {trk)k a b - 87rV a 0V fc( /.. 

Here R is the scalar curvature of g ab and R ab its Ricci tensor. It follows from 
the Einstein-scalar field equations as a consequence of the Bianchi identity that <f> 
satisfies the wave equation W g^X} ' a V 'p<j> = 0. This has the 3+1 form: 

-<9 t 2 0+(trfc)<9 t + A</> = O. 

The constraints and evolution equations are together equivalent to the full Einstein- 
scalar field equations. 

As mentioned in the introduction, generic non-oscillatory singularities are ex- 
pected to have the feature that spatial derivatives become unimportant near the 
singularity. We will make an ansatz for g ab , k ab , <fi which has leading order terms 
°gab, kab,°4' solving & version of the Einstein evolution equations and the wave 
equation with (most) spatial derivatives cancelled. This is the velocity dominated 
system, with constraint equations 

-°k ab k ab + (tr k) 2 = 8ir(d t °<j>) 2 , 

V Q (°fc a6 ) - e 6 (tr°fe) = -8™> t Ve 6 (V), 

and evolution equations 

d t °gab - -2°k ab , 

d t °k a b = (tr°k)°k a b - 8nd t °(t) 2 S a b . 

The velocity dominated scalar field °4> satisfies the equation 

-9 t 2 (°0) + (tr°A:)9 t V = O. 

The velocity dominated evolution equations are ordinary differential equations, 
while the velocity dominated constraint equations still include partial differential 
equations. The space of solutions of the velocity dominated constraint equations 
can be proved to have the same number of free functions as the standard Einstein 
constraint equations. 
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The solution of the velocity dominated evolution equations is, after an ap- 



propriate choice of time, given by k a b 



-t K a b , where K a b (x) is a time- 



independend symmetric matrix, which is assumed to have positive eigenvalues 
{p a (x)}, the Kasner exponents. The velocity dominated scalar field °4> is of the 
form °4> = Alog(t) + B, where A, B are time-independent functions. 

There is an important technical problem when the matrix K a b has double eigen- 
value at some point x. In this case, it is impossible in general to choose a smooth 
frame {e a } which diagonalizes K a b . This causes some difficulties which which we 
will gloss over in the discussion below. 

Given the solution to the velocity dominated system, we make the following 
ansatz for the solution of the full 3+1 Einstein scalar field system. We work in terms 
of a frame {e a } which (approximately) diagonalizes K a b , let q a be (approximate) 
|AlO0| , §5] for details. 

= max(s,0), and let ao be a small positive number, to be 
ao- The ansatz we use is 

- °9ab + WVfc. 

2 (°k C b +t- 1+aC ^ C b ), 



eigenvalues of K a b , see 
For s e K, let (s)+ 
chosen. Let a a b = (qb - 



9ab 
k a b 



9ac\ 



where 7 c b , K c b ,ip are o(l). The ordo notation used here, and in Theorem 2.1 below, 
has a precise technical meaning, cf. [ AROC , §4] . It is important to note that g ab is 
not a priori symmetric but this has to be shown as a consequence of the equations. 

Rewriting the Einstein and scalar field evolution equations using the ansatz gives, 
among other equations, the following evolution equations for 77., n r 

td tl a b + a a bl \ + 2k\ + 2 7 Q e (i°fc e b ) - 2(i°r e ) 7 e b = -2t 



fj' n 6' 



7 e K 6> 

td t K a b + a a b K a b - (i°fc a b )(tr K ) = t a «(tr K )K a b +t 2 - a ^( s R a b - %ng ac e c {cp)e b (cp)). 

Here R a b is the Ricci tensor with respect to a symmetrized version of g a b- It can 
be shown that the power of t occurring on the right hand side of the evolution 
equation for 7 a b is positive. 

In order to write the system in first order form, it is necessary to introduce an 
additional field A a bc , defined in terms of the first order spatial derivatives of 7° b . 
Proceeding in the same way for the wave equation results in a first order system 
in the variables ip, u> a , x, where w a , x encodes first order space and time derivatives 
of ip, respectively. Let u = (j a b , ^ a bc i K<1 b> ^1 x) be the unknown in the resulting 
first order system. The Einstein scalar field system then takes the Fuchsian form 

(2.1) td t u + AU = f[t,x,u,u x ], 

where / = o(l), in the sense that / applied to the fields specified by the ansatz 
tends to zero as a power of t, as t tends to zero, and the matrix function A satisfies 
a spectral condition, related to being positive definite, cf. [AROO, §4]. Under an 
analyticity condition, existence and uniqueness for the system ( 2.l[ ) holds by a 
singular version of the Cauchy-Kowalewski theorem, cf. [ AR00| , Theorem 3]. The 
main technical difficulty in proving that the Einstein-scalar field system is in the 
Fuchsian form ( fO| ) lies in the proof that the function / is o(l) in the appropriate 
sense. For the proof that / is o(l), the most important point is to estimate the 
term t 2 ~ a bS R a b occurring in the evolution equation for n a b . This estimate requires 
a detailed analysis of the blowup rate of the connection forms of the metric g a b, 
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and their derivatives. The Einstein-scalar field and Einstein-stiff fluid systems 
can be written in Fuchsian form only if p a > 0, in the four-dimensional case we 
are considering. The above proves existence and uniqueness of solutions up to the 
singularity, for the Einstein-scalar field system with real analytic data. 



Theorem 2.1 (| ARO0| , Theorem 1]). Let S be a real analytic manifold of dimen- 
sion 3, and let {°g ab (t), °k ab (t), </>(£)) ^ e a rea ^ analytic solution of the velocity 
dominated Einstein-scalar field equations on S x (0,oo), such that ttr°fc = — 1 and 
—t°k a b has positive eigenvalues. Then there exists an open neighbourhood U of 
S x {0} in S x [0, oo) and a unique real analytic solution (g a b(t), k a b(t), 4>(t)) of the 
Einstein- scalar field equations on Un(Sx (0, oo)) such that for each compact subset 
K C S there are positive real numbers £, (3, a a b , for which the following estimates 
hold uniformly on K : 



5\ + o(t a ^), 



°g ac gcb 
k\ = °k a b + o(t- 1+aa "), 

d t ^ = d t ^ + o(t- 1 +f ) ), 

Q g ac e f (g cb ) = o(t aa >-<), 

e a (<t>) = e a (°<P) + o(t0-<). 



The positivity condition on the eigenvalues together with the velocity dominated 
Hamiltonian constraint imply that the function A occurring in °</> has the property 
that A 2 is strictly positive in the velocity dominated solution. Thus vacuum solu- 
tions are ruled out by the hypotheses of this theorem. See [AR00, Theorem 2] for 
an analogous result in the stiff fluid case. See also |Rcn00a| for a discussion of the 
case of Einstein equations coupled to a nonlinear scalar field. It is an interesting 
problem to generalize the results discussed here to C°° data, as was done for the 
Gowdy problem in [ RcnOOt ] . 
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